Energy relaxation in the spin-polarized disordered electron liquid. 
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The energy relaxation in the spin-polarized disordered electron systems is studied in the diffusive 
regime. We derived the quantum kinetic equation in which the kernel of electron-electron colli- 
sion integral explicitly depends on the electron magnetization. As the consequence, the inelastic 
scattering rate is found to have non-monotonic dependence on the spin polarization of the electron 
system. 
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Solid state mesoscopic electronic systems provide an 
outstanding microlaboratory for various experiments on 
low-temperature physics. In particular, it allows to study 
the fundamental relaxation mechanisms in solids that is 
important for finding a way to create advanced cryogenic 
devices. Last decade large progress have been achieved 
in fabrication of devices in which electron distributions 
can be controlled and manipulated. For example, they 
are electronic refrigerators, thermometers, radiation de- 
tectors, and distribution controlled transistors The 
main building block of the devices is a diffusive normal 
metal or a heavily doped semiconductor wire connected 
to massive electrodes acting as reservoirs. One of the 
basic goals of the manipulation with electron distribu- 
tions is to cool electrons in the diffusive metal or the 
wire much below the lattice temperature. In general, 
electrons interact among themselves; also they are cou- 
pled to phonons and to the electromagnetic environment. 
Qualitative understanding of the performance of the de- 
vices is based on the ability to solve the quantum ki- 
netic equation for the electron (spin) density matrix / in 
which interactions are taken into account via scattering 
integrals [2J, l3j ■ At low temperatures (T) which are typi- 
cal for experiments electron-electron interaction provides 
usually the strongest mechanism for energy relaxation. It 
can be characterized by the out-scattering rate (l/r ou t) 
that appears in the kinetic equation formalism and has 
the meaning of the inelastic electron-electron collisions 
frequency [4|. With its help one can estimate time and 
length scales at which a non-equilibrium electron distri- 
bution function, e.g., after the cooling cycle, can be ap- 
proximated by the Fermi-Dirac distribution with some 
effective electron temperature 

Usually, a magnetic field is used as one of conve- 
nient tools for manipulation with the electron distribu- 
tion function . In a weak magnetic field (the Zeeman 
energy is much smaller than the Thouless energy) the 
collision integral /[/] acquires magnetic field dependence 
only due to the shift of electron energy in the distribu- 
tion function due to the Zeeman splitting [7]. However, in 
typical experiments [J, [1] the Thouless energy is compa- 
rable with the temperature: T ~ 0.1 K, and, therefore, 



magnetic field above 0.1T should be treated as strong. 
In this regime the magnetic field not only shifts elec- 
tron energy in the distribution function but changes the 
frequency dependence of the scattering probability such 
that a kernel of the collision integral becomes magnetic 
field dependent. This effect should be especially impor- 
tant for the disordered electron systems due to diffusive 
propagation of electron- hole excitations. 

In the Letter we consider the energy relaxation due to 
electron-electron interaction in disordered electron sys- 
tems in the presence of a large Zeeman splitting 2|m|/z/, 
where v denotes the thermodynamic density of states per 
one spin projection. The average spin density m of the 
electron system caused either by the applied magnetic 
field or by an appropriate boundary conditions. We de- 
rive the quantum kinetic equation in which the kernel 
of the inelastic collision integral explicitly depends on m 
(cf. Eqs. l[T8] ) -([T9] ) ). Wc find that in the out-scattering 
rate which can be written as the sum of contributions 
from the singlet (I/t^) and triplet (l/ T out) channels: 
l/T ou t = 1/t-W + 1/7-^(0) + 2/r^ t (|m|), contributions 
from the triplet channel with non-zero spin projection 
on the m direction are strongly affected by the Zeeman 
splitting. For example, at T = and in the case of d = 3 
dimensions we obtain (we use units H = c = ks = 1) 
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Here, 7 = -Fj{/(1 + Fq ) where F£{F£) stands for the 
standard Fermi liquid interaction parameter in the triplet 
(singlet) channel, D denotes the diffusion coefficient, and 
e a = e— ojm|/i/ where e is an electron energy with respect 



to the Fermi energy Ep. 
following asymptotics 
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Therefore, at zero temperature the 1/t 0ui 
non-monotonic function of |m| at fixed quasiparticle en- 
ergy e CT . It has the maximum when 2|m|/^ ~ e CT , and is 
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FIG. 1: (Color online) The color density plot of 
T-<£t(0)/Tout(l m l) as a function of Fq and 2|m|/(i/e tr ). 



strongly suppressed at larger |m| (see Fig. |T|). This non- 
monotonic dependence of 1/T„*t(|m|) on |m| is similar to 
the non-monotonic behavior with magnetic field for the 
scattering rate of electron on a magnetic impurity [jj [Toj ] . 

The most general method to describe a non- 
equilibrium low-energy dynamics in disordered inter- 
acting electron systems is the Keldysh non-linear a- 
model fjjj]. In the presence of non-zero Fq as well as 
m, one can follow the same derivation of the quantum 
kinetic equation as given in Ref. 



Ill, OJ, [13( except for 



complications that arise from non-commutativity of dif- 
ferent components of plasmon and gauge fields. Below 
we sketch the derivation of the quantum kinetic equation 
briefly and highlight the points where it differs from one 



of Refs. 11, 12, 1 



We write the grand partition function of the inter- 
acting electrons in a random potential in the coher- 
ent state basis: Z — J DtpDip exp{iS[ip,ip)}, where 
Sfy, V] = j* c dt { J x sp(^d t V>) - H$, ip}}. Here, C is the 
Keldysh contour, the symbol sp denotes trace over spin 
degrees of freedom and H = Ho + i?i n t ■ The one-particle 
Hamiltonian Ho involves the par ameters with standard 
Fermi liquid renormalizations [14| . The interacting part 
Hmt = J dvdv' {ip r r s (r - r')p T , + 2m r r t (r - r')m r '} 
where T t (q) = Fg/(2v) and T s (q) = V (q) + Fg/(2v) 
contains the long-range part of the Coulomb interac- 
tion Vo 14, [TEI]. The charge and spin density oper- 
ators are given as p r — ^ ia '4>a(rt)ip t j('Ft) and m r = 
| J2 a Vv(rt)S(T(T'Vv(rt), respectively, where the Pauli- 
matrices s a , a — 0, 3 act in the spin space. 

To derive the non-linear a-model we perform stan- 
dard steps [ll|: i) we average Z over the Gaussian, S- 
correlated disorder and then, introduce the Q-matrix; 
ii) next, we decouple four-fermion interaction terms in 
-ffint by the Hubard-Stratonovich transformation using 



0, 3 ; iii) then, we perform the 
Keldysh rotation [l6| and iv) finally, we integrate out 
fermion degrees of freedom. Hence, we obtain 



iS[Q,&] = Trln 



id t -£+-Q + (Q<*-<i><*)s a * ij 



- dtdv[e T T- l a z Q]-^-Tr(Q 2 ), (3) 



where symbol Tr denotes the trace over the spin and 
Keldysh spaces combined with time and space integra- 
tions, r is the elastic scattering time. The Pauli matrices 
a z , 7x = (To and 72 = cr x operate in the Keldysh space, 
£ = p 2 /2m e — Ep where m e is the electron mass, and f is 
the diagonal matrix in the 4x4 9-space: T — F/(2v) = 
diag{r s , r t , r t , r t }. We assume that there are an elec- 
tric potential ip and a static magnetic field H . Then, the 
classical components 4»f =0 — &p and </>" >0 = —gpsHa 
where g and ps denote g- factor and the Bohr magneton, 
respectively; the quantum components 4>2 = 0. The low- 
energy description is valid under the following conditions: 
< 1 and Tt « 1 « E f t, see also Figf2] In addi- 
tion, we assume T/Ep < reH /m e <C 1 and, therefore, 
ignore the Cooper channel and orbital effects. 

The O-field has a nonzero (zero) average for the clas- 
sical (quantum) components, [(. . .) = /(. . .) exp(iS')]: 



<er°> 



<e?>°) = -^m Q , (9?> = (4) 



where p is the average charge density. From here on- 
wards, we omit Vb(q) for a sake of simplicity: it can be 
always restored by taking the so-called "unitary" limit, 
Fq — > 00 Next, to improve the convergence of expan- 
sion of the logarithm in Eq. ([3lj , we perform the gauge 
rotation of the Q-matrix 



Q t , v (r) = U{t,T)Q ttt ,{v)U-\t',r), Q=(^ K A j 

U(t,v) = exp{i{k)} exp{iSk} , k = kf^ l s a . (5) 

In general, the gauge k and plasmon 8 fields can 
be separated into slow, (k) and (0), and fast, Sk and 
SQ, contributions as compared to the Thouless en- 
ergy (see Figf2]). Physically, the fast components de- 
scribe charge and spin fluctuations in the electron sys- 
tem. Under assumptions that Ep 1/r 3> T and 
r s = e 2 /vp < 1 where v F is the Fermi velocity, we ex- 
pand the action ([3]) around the standard saddle-point 
Q = A with Z = 0, R = —A = 1 and arbitrary K. 
In the equilibrium, the Wigner transform K e (r,t) — 
J dt'K t+tl/2 , t -t'/2('r)(exp(i€t') is equal to K^ q (y,t) = 
2tanh[(e — <p — dt(kf)s a )/2T}. As usual, we restrict our- 
selves to the second order in Sk, SQ and gradients of Q 
within the low-energy manifold Q 2 = 1 . 
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FIG. 2: "Large" and "small" energy scales of the problem. 
Here S is the mean level spacing and L is characteristic sys- 
tem length. The slow in time fields, e.g. (8 a ), have the 
Fourier transforms which are localized within the energy in- 
terval (cT,a>,nin)' The Fourier transforms of the fast in time 
fields, e.g. SQ a that describe the density and magnetization 
fluctuations, are inside (u m i n ,w ma!t ). Similar diagram can be 
drawn for the definition of slow and fast spatial scales. Inte- 
grating out fast degrees of freedom in the partition function, 
Z, we obtain the kinetic equation for the Wigner transform 
of the one-particle density matrix which is slow in space and 
time. 



Then, we find 

iS = i J dt J dr[e T f- 1 CT x 6 + vb T a x b\ 

-^-[D Tr{d r Q) 2 + U Tr(id t + b)Q] . (6) 

Here, d r Q = VQ + i[g, Q]_, b = U-\-4> + Q)U + 
U~ 1 [id t , U]_ and g = -eA/c + U^ 1 [p, J7]_ where A 
stands for an external vector potential. For latter conve- 
nience, we introduce the notations: b = (fy + S^'b+S^b, 
where (b) — (0) — c/> — dt (k) is of the zeroth order in 5k 
and 50, S^b is the first order term, and S^b is of the 
second order. 

The particle and spin densities can be found from 
Eq. © as 
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tv(a x s a (Q t , t )) + 2is(b)«, (7) 



where the trace is over the Keldysh and spin spaces. 
Given Eqs. (|4]) and (0, it is trivial to derive 



P = 



■Kvti{a x s (Q)) 
2(1 + Fg) 



2v 



-[tp + dtikl)}, (8) 



irvix(& x s a {Q)) v[gn B H a /2-d t {k^)} 
m = — 1 : — ■ (9) 



4(1 + F%) 



1 



We emphasize that the interaction renormalizations of 
charge and spin densities in Eqs.fJS])-© are in agreement 
with the Fermi- liquid theory 14 1. 

Although the theory ([6]) encodes all low-energy dynam- 
ics of the electron system, for accurate derivation of the 
quantum kinetic equation it is enough to consider only 
a saddle-point configuration Q for a given configuration 
of plasmon and gauge fields. The saddle point (Usadel) 
equation is as follows 



Dd r {Qd r Q)-[dt-tb,Q}- =0. 



(10) 



A general solution of the Usadel equation can be written 
as Q = (Q) + 5Q where (...) denotes the average over 



SO fluctuations from here onwards. The term SQ involves 
fluctuations <50 and 8k of the plasmon and gauge fields 
that are controlled by the small parameter (Epr)^ 1 <C 1. 
Given Eq. (fT0|) . it is sufficient to find SQ to the second 
order in 50 and 5k. This procedure can be significantly 
simplified if we demand that Z-component of 5Q does 
not contain linear in 50, 5k terms. It is so if 



{ (1) tJ-J)V{ (1) g 2 = 0. 



(11) 



Due to the nonlinear condition Q 2 = 1, this relation au- 
tomatically ensures that retarded (R) and advanced (A) 
components of 5Q have the same smallness as the Z- 
component. We find the relation between classical com- 
ponents of 50 and 5k by requiring the Keldysh (K) com- 
ponent of the Usadel equation with {Q} substituted for 
Q to vanish in the linear order in 6k and 50 at t — > t': 
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■ 5b)? (w) = -2w[5- 1 J B a) ]^£>V5gf (w) . (12) 
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(13) 
(14) 



with e± = e ± lo/2. In the absence of magnetic field, and 
in the equilibrium, we have B^f = S a/3 coth(w/2T) and 
= aj5 af3 . Then, Eqs. (JT]) and (UJ) reduce to the 



corresponding conditions of Refs. [llj, [12|, |13| . In general, 
we find = ujI + A, where 

de 
2^ 



X 



a/3 



we** 



-<*«>} 



(15) 



Derivation of the quantum kinetic equation becomes 
less cumbersome in the (fc)-gauge: (b) = 6. Then, it 
implies tfWg = VSk, S^b = 50 - d t 5k and S^b = 
i[50 — TjdtSk, 5k]-. We mention that in the (fc)-gauge A Q ^ 
intimately related with the average spin density: \ a P = 
2ie a p 1 m y /v. It is the presence of non-zero A Q/3 , S^b and 
the matrix structure of B^j 3 that strongly complicates the 
derivation of the quantum kinetic equation for m^O. 

Substituting Q into the action we find after ex- 
pansion to the second order in 50: 

iS[50] = -iiruTi(6^b 2 (K))+i J cttdr[c50 T f- 1 a x 50 

+v5^b T a x 5^b] + ^ Tr[5«g, A] 2 .. (16) 

Given Eq. (fl6|) . it is easy to find the 2-point correlation 
function V°f(rt,r't') = i4vD(6^sf (ri) $( x )gf (r't')) as 



Vn(q,uj) 



Dq 2 
Dq 2 — iuj 



Vx2(q,u) = -2iuJVii(q,w)[B u - i(Dq 2 - + F^ 1 ) 
xS-'B^ + i(Dq 2 + iu)BJz X {1 + F- x )]V 22 {q,u), 
V 22 (q,uj) = [V xl {q,w)}\ V 21 (q,iu) = 0. (17) 
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Substituting Q — (Q)+SQ into Eq. (fT0|) . then, expanding 
its left hand side to the second order in <50, and, finally, 
averaging the result over the <50-fluctuations with the 
help of Eqs. (fTT|) . (fl~2|) . and (jTTJ) , we obtain the quantum 
kinetic equation for the one-particle spin density matrix 
f=[2s°-(K e (v,t))]/A: 



DAf - d T f 



, E + V(m ' s) )Dd e Wf = I[f], (18) 



where E stands for the electric field, ft is convenient to 
choose z axis along m, then / becomes diagonal and the 
collision integral acquires the following form: 

i[U] = —l ^{[pWKo) + pW(^o)R„(^) 

V J Z7T I 

(19) 



+2P< ff >(w,|m|(7).V a (e,w)}. 



Here, 

Ja ><T >(e,w) 



d£_ 
2^ 



fe 'aft 



(1 - fe' ,ct)(1 - fe-w,(T')fe'_,(7'fe,<7 



(20) 



pW( W ,H)=x; 



,(i-/ £ ,_ )CT/ )(i-/ ei(T ) 

(F$Dq 2 \Dq 2 - iuj^ 1 ) 2 



\D(1 + F«)q 2 - i{uj + 



Due to commutativity of the total spin with the Hamil- 
tonian, /[/] does not lead to "spin- flip" processes; the 
spin density evolves according to the equation: <9 r m = 
DV 2 m. It is worthwhile mentioning, that in the equilib- 
rium / £ (T becomes the Fermi-Dirac distribution function 
/f(co-) and ^[/<r] vanishes identically. 

The out-scattering rate can be found from I[f a ] by its 
variation over / £ (T at the equilibrium 0, Then, from 
Eqs. m we find l/r out = l/r^ ) t + l/r^ ) t (0)+2/r^ ) t (|m|) 
where 1 /r^ t is given by the standard expression 0, 0| 
and 



_(*)/ 

'out \ 



J 2-kv 

to [coth £p + tanh 



(21) 

with Y(e a ,uj,T) = to [coth jfc + tanh . The 
p( a )((jj + 2|m|er/^, |m|) is determined by the screened 
electron-electron interaction and by the probability for 
electron to diffuse 0]; in agreement with Ref. [l||, they 
involve diffusion propagators for spin excitations in the 
presence of m. Evaluating integrals over momentum and 
frequencies in Eqs. (|2Tj) for T = and d = 3 we obtain 
the result ([1]). As known very well in lower dimensions 
and non-zero temperature one should evaluate the inte- 
grals in Eq. (|2"T|) self-consistently [HI ; we present detailed 
results for T > and d = 1, 2, 3 elsewhere [l9j |. 



In summary, we derived the quantum kinetic equa- 
tion that describes the energy relaxation due to electron- 
electron interaction in disordered electron systems in the 
presence of non-zero spin polarization. We found that 
the T = rate of electron-electron collisions is non- 
monotonic function of the electron magnetization. It can 
be used for decoupling electron degrees of freedom from 
the environment. 
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